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0-( Abstract 

1) : 

f^, ' We consider the Lagrangian of gravity covariantly amended by the mass 

and polynomial interaction terms with arbitrary coefficients, and reinvestigate 
^^ I the consistency of such a theory in the decoupling limit, up to the fifth order in 

^f^ ■ the nonlinearities. We calculate explicitly the self-interactions of the helicity-0 

'sj" , mode, as well as the nonlinear mixing between the helicity-0 and -2 modes. 

We show that ghost-like pathologies in these interactions disappear for special 
choices of the polynomial interactions, and argue that this result remains true 
f— ^ ■ to all orders in the decoupling limit. Moreover, we show that the linear, and 

^^ I some of the nonlinear mixing terms between the helicity-0 and -2 modes can 

be absorbed by a local change of variables, which then naturally generates 
the cubic, quartic, and quintic Galileon interactions, introduced in a different 
k> , context. We also point out that the mixing between the helicity-0 and 2 

^ I modes can be at most quartic in the decoupling limit. Finally, we discuss the 

implications of our findings for the consistency of the effective field theory 
away from the decoupling limit, and for the Boulware-Deser problem. 



o 



c^ 



1 Introduction and summary 

In this work we study the covariant polynomial potential of a relativistic and sym- 
metric rank-2 tensor field living in four-dimensional fiat space-time. 

We start with the mass term in the potential. Poincare symmetry in four di- 
mensions imposes that any massive spin-2 state has to have five physical degrees of 
freedom - namely, the helicity-±2, helicity-±l, and helicity-0 modes. The quadratic 
potential that describes these degrees of freedom is that of Fierz and Pauli (FP), 
[1]. The latter is known to be the unique ghost-free and tachyon-free mass term for 
the spin-2 state [2]. 

No matter how small the graviton mass is in the FP theory, the helicity-0 state 
couples to the trace of the matter stress-tensor with the same strength as the helicity- 
2 does [3]. This discontinuity would rule out, on simple observational grounds, the 
FP mass term for gravity. 

As argued first by Vainshtein, the discontinuity problem can be cured by the 
nonlinear interactions which would become comparable to the linear terms already 
for very weak fields [4J. Then, the non-linearities could give rise to the screening of 
the helicity-0 mode at observable scales, rendering the theory compatible with the 
known empirical data |H |5] . 

However, the very same non-linearities that cure the discontinuity problem typ- 
ically give rise to a ghost in massive gravity, [6]. This ghost, sometimes referred to 
as the Boulware-Deser (BD) mode, emerges as a sixth degree of freedom, that is 
infinitely heavy on a fiat background, but becomes light on any reasonable nontriv- 
ial background {e.g., on a cosmological background [7], or on the weak background 
of a lump of static matter [HI [U [lO]). It is straightforward to see this ghost in the 
so-called decoupling limit [8], in which the dynamics of the helicity-0 mode can be 
made manifest. Then, the sixth degree of freedom ends up being related to the 
nonlinear interactions of the helicity-0 mode [Bl [U llOjH 

The obvious question to ask is then whether there exists a nonlinear model that 
exhibits the Vainshtein mechanism, but without the ghost mode. This question was 
raised in Ref. [8], and studied in detail in Ref. ^. The latter work argued that at the 
cubic order the ghost can be avoided by tuning the coefficients of the quadratic and 
cubic order terms. Recently, the cubic terms were calculated in a nonlinear massive 
spin-2 theory of Refs. [HI [15], where it was shown that the necessary tuning is in 
fact automatic in this model, and the theory is ghost-free to that order p!6]! 

In the present work we focus instead on addressing this question at higher orders, 
and in a model-independent framework. We therefore allow for arbitrary nonlinear- 
ities in the potential up to the quintic order, but restrict ourselves to considerations 
in the decoupling limit only. 



^Notice also that the discontinuity is absent when a smah cosmological constant is included 
before sending the mass of the graviton to zero [11] [12]. Doing so in de Sitter space, however, one 
passes through the parameter region where helicity-0 becomes a ghost [T^ , while the anti de Sitter 
case is ghost- free [HJ [12] . 



Our result clashes with one of the conclusions of Ref. [9] which states that the 
quartic interactions in the decoupling limit ineradicably lead to a ghost. Regretfully, 
the decoupling limit Lagrangian obtained in Ref. [H] is not reparametrization invari- 
ant neither at the cubic nor quartic order, and gives a tensor equation that does not 
satisfy the Bianchi identity. The ghost found in the decoupling limit of Ref. ^ is an 
artifact of these properties. Hence, we re-investigate this issue in the present work. 
We find a decoupling limit Lagrangian that is similar to that of Ref. P], but differs 
from it in detail, by coefficients of various tensorial structures. In particular, due to 
those coefficients, our Lagrangian is reparametrization invariant, and naturally leads 
to a tensor equation for which the Bianchi identity is automatically satisfied (as it 
should be since the helicity-2 mode only mixes linearly in the decoupling limit). 
Then, not surprisingly, we arrive to a different conclusion, that the quartic theory 
is also ghost-free in the decoupling limit. Moreover, we go on one step further and 
investigate the quintic-order theory, which we also show is ghost-free in the decou- 
pling limit. This also allows us to understand the structure of the interactions to 
all orders and to argue that the decoupling limit can be at most quintic order in 
interactions (or quartic in the mixing between the helicity-0 and 2 modes) in the 
ghost-free theory. 

Finally, as a corollary, we find that the decoupling limit of the most general 
consistent theory of massive gravity gives rise to the quadratic, cubic, quartic and 
quintic Galileon kinetic interactions introduced in Ref. [17] in a different context 
(namely, as a generalization of the special cubic term appearing in the decoupling 
limit of DGP [H] found in Ref. [19]). The Galileon interactions share the important 
properties of (i) being local, (ii) preserving the shift and galilean symmetry in the 
field space of the helicity-0 mode (in particular, in the kinetic and self-interaction 
terms but not in interactions with matter), (iii) giving rise to equations of motion 
with a well-defined Cauchy problem. Since then, the Galileons have developed 
their own independent and interesting life (see, e.g., [201 EI])- We show here that 
the Galileons naturally arise in the decoupling limit of a general theory of massive 
gravity. This also helps to prove that upon appropriate choices of the coefficients 
in the potential, the decoupling limit of massive gravity is stable, at least up to the 
quintic order in interactions. 

We continue this section with a discussion and summary of our main results in 
more technical terms, before turning to the detailed calculations in the subsequent 
sections. 

In analogy with a massive non-Abelian (Higgs-less) spin-1 [22], the dynamics 
of the helicity-0 mode, tt, can be extracted in a generic theory of gravity with a 
nonlinear potential by taking the decoupling limit [8] 

m->0, Mpi ^ oo, keeping As = (m^Mpi)^/^ fixed. (1) 

Following ^, in a generic case of the nonlinear potential, the corresponding La- 



grangian for the helicity-0 mode reads schematically as follows: 

.,.!.□., (^. (2) 

The cubic interaction with six derivatives gives rise to a ghost on locally nontrivial 
asymptotically-flat backgrounds {e.g. on the background of a local lump of matter). 
This could be seen by observing that for vr = tt^' + 6tt, with tt'^'- denoting the weak 
field of a local source, and 611 its fluctuation, the cubic term in (^ could generate a 
four- derivative quadratic term for the fluctuations. This leads to a ghost, which is 
infinitely heavy on Minkowski space-time, but becomes light enough to be disruptive 
once a reasonable local background is considered, see Refs. [H HI [TO] . 

To avoid pathologies such as in ([2]), the Fierz-Pauli combination in the graviton 
potential should be pursued further by tuning the coefficients of various higher order 
terms. This leads to a cancelation of all the terms for vr that are suppressed by the 
scales A5, A4 = (m^Mpi)^/^, A11/3 = (m®M|i)^/" etc. ..for any scale A < A3 = 
{m'^MpiY^^, such that only the terms suppressed by the scale A3 survive. Then, A3 
is kept fixed in the decoupling limit, and the surviving terms (in addition to the 
linearized Einstein-Hilbert term) read as follows: 

A^ = h^' (41^ + i|4l^ + Jl^lS) . (3) 

Here, h^u denotes the canonically normalized (rescaled by Mpi) tensor field pertur- 
bation, while XlJ,XlJ, and X^J are respectively, linear, quadratic and cubic in vr. 
Importantly, they are all transverse (for instance, X^J oc tj^J^tx — d^d^n). Not only 
do these interactions automatically satisfy the Bianchi identity, as they should to 
preserve diffeomorphism invariance, but they are also at most second order in time 
derivative. Hence, the interactions ([3]) are linear in the helicity-2 mode, and unlike 
the previous results in the literature, present perfectly consistent terms, at least up 
to the quintic order. 

Furthermore, some of the terms in (|3]) can be absorbed by a local field redefi- 
nition. For instance, the quadratic term, h^'^X^J, can be absorbed by a conformal 
transformation /i^,^ — )■ /i^,^ -|- r/^jyVr. This shift, besides removing the above mixing, 
generates terms of the form vrX^^^ and 7rX*^^\ which coincide, up to a total deriva- 
tive, with the cubic and quartic Galileon terms [IT]. Further diagonalization of the 
cubic mixing term, h^'^Xf^u , also generates the quintic Galileon, hence exhausting 
all the possible terms that can arise in the Galileon family at arbitrary order. 

Moreover, we also point out that if the decoupling limit happens to pick up the 
scale A3 (as opposed to another smaller scale such as A5, A4, etc. . . ), the mixing 
between the helicity-0 and -2 modes must stop at the quartic order. Therefore, 
for appropriate choices of the interaction coefficients, the decoupling limit at this 
order is exact! It is the subsequent diagonalization of the nonlinear terms in the 
Lagrangian that generates the quintic Galileon. 



Finally, the absence of a ghost in the decoupling limit does not prove the stability 
of the full theory away from the limit and the Boulware-Deser ghost is still expected 
to be present in general. However, it at least shows that one has a well-defined 
and consistent effective field theory below the scale A3. Above this scale, the full 
theory has to be specified. We discuss related issues in section 5. Before that, our 
work has a two-fold motivation: (i) To establish a consistent effective field theory 
below A3 (for the full theory to be viable its decoupling limit should be ghost-free 
as a necessary condition), (ii) All the known examples show that the Boulware- 
Deser ghost, if present in the full theory, does also show up in the decoupling limit. 
Therefore, it is encouraging to find no ghosts in this limit. 

The paper is organized as follows: In section 2 we summarize the formalism 
used to study the decoupling limit of massive gravity with a general potential. We 
then explicitly compute the decoupling limit Lagrangian to the quartic and quintic 
orders in section 3. We work with a generic nonlinear completion of the FP gravity 
for which the scale A| = Mpim'^ is fixed. We argue that the vr mode does not 
decouple from the tensor mode, but that the interactions are free of any ghost-like 
pathologies. In section 4 we give a general framework for computing the Lagrangian 
in the decoupling limit, and argue that in theories which are consistent with the fixed 
scale A3, at most the quartic order mixing term can be obtained, all the higher order 
mixing terms being zero. Moreover, we show in section 5 that upon an appropriate 
change of variables we recover the standard Galileon interactions. Section 6 contains 
some discussions on open issues and future directions addressing the consistency of 
massive gravity away from the decoupling limit. 



2 Formalism 

2.1 Gauge invariant potential for gravity 

Below we consider in detail the decoupling limit of a general Lagrangian of a massive 
spin-2 field endowed with a potential on Minkowski space-time. We use the technique 
developed in Ref . [8] . The covariant Lagrangian with the potential reads as follows: 

C = M^,y^R-^^y^{U2{9,H) + U,{g,H) + U,{g,H) + U,{g,H) ■ ■ ■) , (4) 

where f/j denotes the interaction term at i*^ order in H^^, 

U2{g,H) = Hl-H\ (5) 

U^{g,H) = c,Hl + C2HHl, + csH\ (6) 

U,{g,H) = d,Hl^ + d^HHl + d,Hl^Hlp + d,H^Hl^ + d,H\ (7) 

U,{g,H) = hHl + hHHl + hH^Hl + UHl^Hl, (8) 



Here the index contractions are performed using the inverse metric, so that H = 
g'^'^H^y, H'^^ = g^^'^g°'^H^aH^i3, etc. . . . The coefficients Cj, di and /, are a priori 
arbitrary, but will be determined by demanding that no ghosts are present at least 
up to the quintic order in the decoupling limit. 

Finally, the tensor Hf^iy is related to the metric tensor as follows: 



Mpi 

= H^, + riabd^rdu(t^\ (9) 

where a, 6 = 0, 1, 2, 3, rjab = diag(— 1, 1, 1, 1), and H^y is a covariant tensor as long as 
the four fields 0" transform as scalars under a change of coordinates. Furthermore, 
expressing (()"■ in terms of the coordinates x", and the field vr" as cl)"" = (x° — vr") 5^, 
we obtain 

H^,u = -TJ^ + df,TTy + d^TCf, - r]ai3df,TT°'dy'!T'^ . (10) 

Mpi 

In flTOj) . and in what follows, we adopt the convention that the indices on vr^ are raised 
and lowered with respect to the Minkowski metric rj^,^. Crucially, the expression for 
the tensor ff^^ in ( ITOl) differs by a minus sign in front of the last term from the 
analogous expression in eq. (5) used in Ref. [9j. To emphasize the importance of 
this sign, we derive in Appendix A| the decoupling limit using the opposite sign in 



(TTOj) . recover the results of Ref. [9j, and show that the Bianchi identity is then not 
automatically satisfied, since the reparametrization invariance is not retained in the 
resulting Lagrangian. 

From (jll) it is not immediately clear what is the scale of the effective field the- 
ory represented by this Lagrangian, i.e., what is the energy /momentum scale by 
which the higher polynomial interactions would be suppressed as compared with 
the leading ones. This will become clear by studying the decoupling limit of the 
theory. 

In what follows, we focus on the helicity-2 and helicity-0 modes, but ignore the 
vector mode. The latter enters only quadratically in the decoupling limit (since the 
vector does not couple to a conserved stress-tensor in the linearized order), and can 
be set to zero self-consistently. This does not prove that the vector sector is ghost- 
free, however, the findings of Ref. [16] that the cubic nonlinear it ies for the vector 
are completely harmless due to the U{1) gauge invariance of the resulting terms, 
suggest that the vector sector is not going to reintroduce the BD ghost. Therefore, 
we use the substitution: vr^ = c^aVr/Ag, so that 

Illy rO / \ 

where we use the same notation as in [_9j, 11^,^ = d^dyj: and 11^^ = rj'^^Ii^aJ^py. 
Moreover, in what follows the square brackets [. . .] will represent the trace of a tensor 
contracted using the Minkowski metric, e.g. [11^] = II^^II^j, and [11]^ = II^IIJ^. 



2.2 Decoupling scale 

As mentioned in the introduction, the interactions f/2 and f/3 typically lead to terms 
of the form ((9^7r)^/(Mpim^), and the decoupling limit should be taken keeping the 
scale A5 = Mpim^ fixed, while Mp\ — )■ 00 and m — )■ 0. However we will show in 
what follows (see also [16]) that for some special values of the coefficients q, such 
interactions cancel (up to a total derivative), generalizing the FP term to the cubic 
order. This procedure can be extended further to an arbitrary order: 
At a given order the leading contributions are of the form 

{ddnr .... 

then, one chooses the interactions Un{H) ~ H^ so that the above terms combine into 
a total derivative. At each order, there exists a unique total derivative combination 
'^der t^^* '^^^ t)e written as follows: 



m=l ^ '' 

with £|jgj. = 1 and £jjj. = [11]. Up to the quartic order, the total derivatives are 

4el = [n]^-[n^ (14) 

4'el = [nf-3[n][tf] + 2[n3], (15) 

£W = [n]4-6[n2][n]2 + 8[n3][n] + 3[n2]2-6[n^]. (16) 

Moreover, at higher orders these total derivatives vanish identically, i.e. Juf^^ = 0, 
for any n > 5. By ensuring that all the leading terms flT2|) take the form of a 
total derivative flT3|) . all the interactions that arise at an energy scale lower than 
A3 disappear. Keeping this in mind we will therefore consider below the following 
decoupling limit (firts considered in [19] in the context of the DGP model) 

m -^ 0, Mpi -^ 00, keeping A3 = (m^Mpi)^/^ fixed. (17) 

Note that the procedure of taking the limit in the present case is well defined for 
fields that decay fast enough at spatial infinity. For these we introduce an infrared 
regulator of the theory, say a large sphere of radius L ^ 1/m, and take the radius 
to infinity, L — )■ 00, before taking the limit (IT7|) . This hierarchy of scales enables us 
to put all the surface terms to zero before taking the decoupling limit. 

Furthermore, as it should be becoming clear from the above discussions, the 
scale A3 will end up being the effective field theory scale. The higher interaction 
terms, both written or implied in (jl]), will be subdominant to the leading ones for 
energy /momentum scales below A3. 



3 Decoupling limit of massive gravity 

3.1 Cubic order 

We now explicitly compute the decoupling limit for the interactions considered in 
dMl]), and thus generalize the Fierz-Pauli term to higher orders. In terms of the 
"Einstein operator" £ defined for an arbitrary symmetric field Z^^, as 

S^^^Z^p = -1 {UZ^, - d^d^Z: - d^do^Zl + d^d^Zl - i]^.UZl + i^^^d^dfiZ'^^) ,(18) 

the decoupling limit Lagrangian of massive gravity up to the cubic order reads as 
follows 

L = -h^^Sf^Kp + h^^X^^iJ (19) 

4 ((8ci - 4)[n=^] + (8c2 + 4)[n][n2] + ^c,[nf) + ^.h'^^xfi^ 



with 



4^5 vv-i -n^^ J > V--. . -yi^^jL^^ j > --.l^^j j ' ^3- ^>. ' 



X« = Mr^.u-n.u, (20) 



and XfiJ quadratic in 11. Using the total derivative combination (fT5|) . the interac- 
tions arising at the scale A5 can be removed by setting 

ci = 2c3 + - and Ca = -Scs - - . (21) 



As a result, we find the following expression for the tensor X^ 



(2) 



xjS = -(6c3 - i){(nj. - [n]n,.) - i ([n^] - [n]^) r^,.} . (22) 

Notice that both XfiJ and Xfi,J are automatically conserved, as they should for the 
reparametrization invariance to be retained and the Bianchi identity to be satisfied. 

Moreover, it is straightforward to check that these cubic interactions bear at 
most two time derivatives, and are therefore free of any ghost-like pathologies. One 
should also check that the lapse (which coincides with hoo in the decouphng limit) 
still propagates a constraint, which is indeed the case here as neither Xqq nor Xqq 
contain any time derivatives. Furthermore, these cubic interactions with the specific 
coefficient C3 = 1/4 have already been discussed in detail in Ref. [T6] . 

We now apply the same formalism to quartic interactions for which ghost-like 
pathologies have been argued to arise inexorably in Ref. [S]. 



d2 -- 


= 8d5 - -{6C3 + I) , 


ds - 


= 3rf5-^(12c3 + l) 


di = 


= -6^5 + ^C3 . 



3.2 Quartic order 

At the quartic order, we find the following interactions in the decoupling limit: 

^^'^ = ^h^'XJ^J + ^{(3ci - 4rfi - i)[n^] + (C2 - 4^3 + i)[nT (23) 

+(2c2 - 4ci2)[n][n3] + (3c3 - 44)[tf][n]2 - 44[n]^} , 

with A4 = {Mpim^Y^^ and XfiJ cubic in 11. Here again the pathological terms 
arising at the scale A4 can be removed by using the total derivative combination 
(irB]) . and by setting ci and C2 as in fl2T]) . as well as 

di = -6ci5 + 4(24c3 + 5), (24) 

lb 

(25) 

(26) 

(27) 

(3) 
Substituting these coefficients in X^^J we obtain the mixing term between the helicity- 

and 2 modes determined by 

xj^J = (C3 + 8d,) {enj, - 6[n]nJ, + sm' - [u'])u,, (28) 

This expression bears two expected but important features: 

(3) 

• It is conserved d'^X^iJ = 0, as it should be for the reparametrization invariance 
to be present and the Bianchi identity to be automatically satisfied. 

• For i,j space-like indices and time-like index: 

(3) 

X!>/ has at most two time derivatives, 

(3) 

Xq^ has at most one time derivative, 
Xqq has no time derivatives. 

These properties ensures that no ghost-like pathology arise at the quartic level in 
the decoupling limit as long as the interactions come in with the generalized FP 
structure set by the coefficients (121]) and (12411271) . 



3.3 Quintic order 

At the fifth order in the decoupling limit, we consider interactions as given in (Q. 
The pathological terms that scale as 

can be canceled with an appropriate choice of the coefficients /i to /g: 

/i = M + 1^3 - 6^5 + 24/7 , /2 = -i - i|c3 + 64 - 30/r , 

/3 = fca - 54 + 20/7 , h = ~TE- 1^3 + 54 - 20/7 , (30) 



/s = ^C3 - 34 + 15/7 , /e = 4 - 10/7 



16 4" 

YgCs - 0U5 -h ±oj7 , /e = 4 - -LUJ7 

As a result, the quintic interactions in vr arrange themselves to form the expression 
for vCdgj., as derived from flT3|) 

££l = 24[nV30[n][nV20[n3]([n]2-[n2]) (31) 

+i5[n][nY - io[n2][n]3 + [n]^ = 0. 

Notice that C\J^ is not simply a total derivative as for the previous orders, but 
instead vanishes identically. This implies in particular that any limiting Lagrangian 
of the form £(") ~ /(n)£|jgj., where / is an analytic function, gives no dangerous vr 
interactions and can be used at higher orders. Beyond the quintic order the degrees 
of freedom in the coefficients to be tuned should therefore increase, and make it 
easier to remove any ghost-like interactions. 

With the above choice of coefficient fl5U]) . the only quintic interaction in the 
decoupling limit then is 

£(^) = Jih'^'Xl'^^ , (32) 

with 

4t^ ~ 24(nj, - nnjj + 124X. - 44tln,. + 4tk. = , (33) 

with £|jg'j. ' given respectively in flTl]) . ( 1TB]) and (TTB]) . The decoupling limit is there- 
fore well behaved up to the quintic order, and the number of free parameters at 
higher orders suggests that one can always make appropriate choices to avoid any 
ghost mode from appearing in the entire decoupling limit. To be certain, one should 
however analyze a fully non- linear theory, such as the one proposed in [HI IS] . 

Motivated by the above obtained results, we set up in the next section a general 
formalism for obtaining the interactions to all orders. 

Before we do so, some important comments are in order. We might of course 
argue that the absence of the ghost up to the quintic order represents no proof of 
the stability of the theory even in the decoupling limit, since the ghost could be 
pushed to the next order in interactions. It is also not a proof of the consistency of 
the full theory, as was discussed in section 1, since the ghost may appear away from 
the decoupling limit. The arguments concerning these two points, respectively, are: 

9 



1. Beyond the quintic order, the number of free coefficients in the interactions 
seems sufficient to ehminate pathological contributions of the form [ddii)"'. 
Furthermore, beyond the quartic order all conserved tensors of the form X/t" ~ 
(557r)" vanish identically, and cannot lead to any ghost-like pathologies in the 

mixing W^Xji, between the helicity-0 and 2 modes. 

2. The ghost may exist in a given order away from the decoupling limit (say at 
the quartic or higher order), but disappear in the decoupling limit. If so, then, 
the ghost should come with a mass greater than A3. Then, the theory would 
be acceptable as an effective theory below the A3 scale. However, at scales 
above A3, one would need to specify an infinite number of terms in the full 
nonlinear theory in order to conclude whether or not the ghost is removed by 
the resummation of these terms. This will be made more precise in the last 
section. 



4 General formulation for an arbitrary order 

All our findings up to the quintic order presented in the previous section can be 
formulated in a unified way, which may also suggest how things could work at 
higher orders. For this, in the iVth order expansion (so far iV < 5), we introduce 
the notations 

M2 rn2 JL 

Ui,{g,H)^ ^-^Y../^U,{g,H), (34) 



4 .., 

where the tensor H^j^y is defined as in section 2. If the A^*'^ order expression for the 
function Uiy{g,H) satisfies 

UN{g,H) = total derivative , (35) 

(where A^ denotes the helicity-1 field) then, the decoupling limit Lagrangian for the 
helicity-0 and -2 interactions, up to a total derivative, takes the form: 

4r = -^^"'^^.t^./^ + h^''Xl!:!\7r) , (36) 

with the conserved tensor X^u '■ 



5h, 



''flV 



(37) 



We have checked that the above Lagrangian gives rise to equations of motion with 
no more than two time derivatives and appropriate constraints for N < 5. It seems 
reasonable to conjecture that this will also be the case for N > 5. Furthermore, 

10 



in four dimensions X^u can only contain a finite number of terms if it is local 
and conserved. It is therefore likely that this formalism leads to a finite number of 
interactions in the decoupling limit. 

At a given order n in the expansion, there should be enough freedom to set the 
polynomial Un{g, H) appropriately, so as to ensure that the leading interactions f lT^ 
enter as a total derivative of the form (iT3l) . or as f{Ii)C^^ for m > 5 and / being 
an arbitrary function of Il^v The resulting leading contribution is then of the form 



£(") = -r^/^^^xW , (3J 



^^3 



where /3 depends on the coefficient c's, rf's, etc. and XliJ ~ H" must be conserved 



(") 



as a straightforward consequence of reparametrization invariance in the decoupling 
limit (since higher interactions in h are then suppressed). At each order ra, there is 
a unique combination of H" 's which is conserved. This combination is of the form 

XfZ^ oc ^^^ . (39) 

In four dimensions however, C\^^ = as pointed out earlier, and the same remains 
true at higher orders. This further implies that there is a limit on the number of 
possible interactions in the decoupling limit: X^Z = for any n > 4. This suggests 
that all theories of massive gravity (with the scale A3) can only have at most quartic 
couplings between the helicity-0 and 2 modes in the decoupling limit. 

5 Massive gravity and the Galileon 



When making the generalized FP choice for the coefficients (1211) . (I24tl27p . and (l3Ql) . 
the higher interactions in the decoupling limit only arise as a coupling between the 
tensor mode and the helicity-0 mode of the form 



>Ci„t = /.-xiiT^ = h^' (x« + j,xj^J + ^x(3)) , 



(40) 



where X*^^^ is given by fl20|) . X^^^ by fl22|) and X^^^ by fl28|) . Moreover, as emphasized 
before, d'^X^l = 0. We proceed further by noticing that 

4'/^ = 4t^a'^ , (41) 

with 



Zi!J = TTT],,, (42) 

Zj^j = (6C3-I)9^7r9,7r. (43) 



11 



We can therefore diagonalize the action up to the cubic order by performing a local 
but nonlinear change of the variable 

h =h + Z'^^^ + —Z^"^^ (AA) 

such that, up to total derivatives, the Lagrangian is 

C = -h^J^^^^^Kp + ^ttDtt + ^ ^^^^^7 ^^ {dirfU-n (45) 

+ X| ^2^^""' " ^^' " ^^'^ + ^^^0 ^^""^^ ^^"'^ ~ ^"^'^ ^ Xf ^''^i'^ 
^ (6c3 - i)(c3 + 84)(5vr)2 ([n]=^ - 3[n][n]^ + 2[n3]) . 



2A1 



In the first line we see appearing the quadratic and cubic Galileon terms, [T7] (the 
usual kinetic term for vr, as well as the interaction present in DGP). In the second 
line we notice the quartic Galileon interaction and finally the quintic, last interaction 
of the Galileon family, appears in the last line. 

By setting C3 = — Sfis we precisely recover the Galileon family of terms up to 
quartic order, and all the remaining couplings with the tensor mode disappear at 
the quintic order. Since there is still a lot of freedom in the coefficients at higher 
orders, it is only natural to expect this result to be maintained to all orders. 

On the other hand, if C3 7^ —8^5, then the last mixing term h^'^XJiJ does not 
seem to be removable via any local field redefinition. This mixing term may be 
crucial to address the issue of superluminality of the massive theory, as the Galileon 
without the mixing terms does exhibit superluminal behavior [17j . 

In a more general case, as soon as the cubic Galileon is present in f l45|) . we are 
also bound to have either the quartic Galileon and no other terms (for C3 = —Sd^), 
or a quartic mixing and the quintic Galileon (for 4(c3 + Sd^) = (6C3 — 1)^ 7^ 0), or 
all of the above terms together. 

If however, the cubic Galileon is absent (for C3 = 1/6), one in general is left with 
the quartic Galileon and the quartic mixing term. 

Finally, notice also that for the specific choice C3 = 1/6 and d^ = —1/48, all 
the interactions at the scale A3 disappear! This may be an example of a theory for 
which the decoupling limit picks up a higher scale A^ > A3, if such a theory exists. 
Alternatively, this may also be a theory in which all the nonlinear terms disappear 
in the decoupling limit. This would suggest that the theory has no strongly coupled 
behavior {i.e. , no Vainshtein mechanism), and would be ruled out observationally. 

6 Outlook 

The previous analysis shows that for appropriate choices of interactions that gener- 
alize the Fierz-Pauli term to higher orders, one can construct a consistent and local 
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theory of massive gravity where no ghost-hke instabihties are present, at least up to 
the quintic order in the decouphng hmit, and positive prospects can be foreseen for 
higher orders. In particular the connection with the Galileon generalization of the 
cubic term appearing in the DGP decoupling limit provides a natural framework for 
studying ghost-free theories of gravity [T71 [21] . 

Furthermore, the decoupling limit considerations of this paper suggest that the 
higher non-linear terms in ([312]) become equally important at the scale A3. Since 
the scale A3 = (Mpim^Y^^ is very low (typically A3 ~ 10~^eV), the effective the- 
ory below A3 can only be used for large scale cosmological studieqj. To extend the 
scope of applicability of massive gravity to shorter length scales, however, one would 
need to go above A3, and, hence, the higher interactions should be taken into ac- 
count. For a viable model, it will therefore be necessary to consider all the higher 
polynomial interactions, Un{g, H), and not only the ones up to the quintic order as 
presented here (even though the decoupling limit may only have a finite number of 
interactions) . 

A theory that provides such a resummation is the model of Refs. [TU [15]. In 
particular, by integrating out the auxiliary dimension in that model, one gets an 
infinite series of interactions of the form ([MH]) and beyond, with certain specific coef- 
ficients. In [TB], it has been checked that the coefficients of the quadratic and cubic 
terms were equal to those used in section 3 for the specific choice C3 = 1/4. Thus, 
in the decoupling limit, the theory is ghost-free up to the cubic order. Furthermore, 
the theory in the cubic order preserves the Hamiltonian constraint even away from 
the decoupling limit [16], and the BD term cancels out in the exact all-order Hamil- 
tonian [H]. Moreover, it was shown in Ref. [I5] that the nonlinear terms giving 
rise to a ghost at a scale A < A3 cancel out in that specific theory. These findings 
constitute an important evidence (but not a proof yet) that the theory of [Tl[ [T5] 
may be consistent, at least classically, to all orders. 

How about other possible theories of massive gravity that would yield the terms 
discussed here with the coefficients still consistent with the absence of the ghost, 
but not coinciding with the ones obtained in [16]? Is there any hope for these 
theories away from the decoupling limit and above the scale A3? Naively, the answer 
seems to be a negative one: As was shown in 0, in the order- by-order expansion, 
and beginning with the quartic order, one cannot avoid higher powers of the lapse 
function in the Hamiltonian, and hence, the emergence of the sixth degree of freedom 
(which typically is a ghost) seems to be unavoidable in massive gravity [9]. 

However, there may be a way to circumvent this problem in the full theory if its 
Hamiltonian, due to a resummation of perturbative terms, ends up having a very 
special dependence on the lapse and shift functions. Here we demonstrate this in a 
toy example, that is motivated by the Hamiltonian of the theory [lU [TJ] discussed 
in [E 



^ Once external classical sources, such as planets, stars, galaxies,.., are present, the energy 
scale of nonlinearities - the Vainshtein scale - depends on the mass/energy of the source and is 
significantly lower |S]. 

13 



Consider the toy Hamiltonian: 

H = N{R'^ + m2/(7)) + Nj {R^ + m^Q'{^)) + m^P{^)^^ , (46) 

where A^, Nj,'yij, and i?", R^ , are the standard ADM variables and functions respec- 
tively [23]; /(7), Qjij) and -P(7) are some functions that modify the GR constraints 
by the mass terms. The shift function Nj is not a Lagrange multiplier, but is alge- 
braically determined, as it should be the case for a massive theory with five degrees 
of freedom. However, the lapse functions also enters in the last term in a way that 
seems to prevent it to be a Lagrange multiplier, and if so, it would give rise to the 
sixth degree of freedom. This is not the case, however: One can introduce a new 
variable rij = Nj /N in terms of which the Hamiltonian reads 

H = A^(i?° + mV(7)) + Nnj{R^ + m^Q^{j)) + Nm^P{^)^— . (47) 

The shift n^-, still has no conjugate momentum, hence 6H/5nj = 0. This determines 
the new shift variable, n^ = —{R^ + m?Q^ {''))) /{m?P {'-))) ^ and yields the following 
Hamiltonian 

Here, the lapse does certainly appear as the Lagrange multiplier. Hence, the BD 
term does not arise, and the theory does not propagate the sixth degree of freedomo. 

On the other hand, a direct perturbative expansion of the last term in fll^ in 
powers oi5N = N—1 with subsequent truncation of this series at any finite nonlinear 
order, necessarily yields higher powers of SN in the Hamiltonian|j. Naively, this 
truncated theory would give rise to the potentially false impression that the lapse is 
not a Lagrange multiplier, and that there is a sixth degree of freedom in the model. 

Noticing that the higher powers of SN at any finite nonlinear order emerge from 
the expansion of the theory ( H6ll is trivial, in this toy model. However, a similar, 
albeit more complicated structure, emerges in the Hamiltonian of the model of 
[T3| [15] (see [H]) and the fact that the terms in the expansion come up from a 
single term in the exact Hamiltonian is not as simple to observe. 



■^In general, it could still be propagating "5.5" modes even if the Hamiltonian constraint is 
maintained. For instance, since the toy model described by (j46]) is not Lorentz invariant for 
general functions /, Qj and P, there may exist non-propagating instantaneous modes in this 
model. For discussions of related issues see, [21]. In contrast, the model of Refs. [TU [TS] is 4D 
Lorentz-invariant and the instantaneous mode in 4D is not expected. For a rigorous proof that 
there are only 5 degrees of freedom, and not "5.5" , however, a detailed study of the algebra of 
the Hamiltonian constraint should be performed. The fact that the decoupling limit gives only 5 
degrees of freedom is an important hint that the full theory is not likely to have the extra "0.5" 
degree of freedom. 

^Note that away from the decoupling limit, and at a nonlinear order, SN is the right variable 
and not hoo = 1 — 7V^ + Nj, which was used before as the lapse in the decoupling limit. 
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Last, but not least, in this work we discussed the classical theory. Generic 
quantum loop corrections are expected to renornialize and detune the coefficients of 
the polynomial terms needed to avoid the ghost. One way to be protected against 
this problem is to have a theory in which the tuned coefficients automatically emerge 
as a consequence of a symmetry that would be respected by the loop corrections. In 
this respect, the recent findings of [16j that the cubic terms with the automatically 
tuned coefficients emerge as an expansion of the theory, which by itself exhibits 
an evidence for a hidden nonlinearly realized symmetry, makes us hopeful for the 
existence of a quantum-mechanically stable effective field theory of massive gravity. 
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Appendix A Decoupling limit with the opposite 

sign in H^^ 

As mentioned in section 12. ![ the expression (1 10 1) for if^,^ differs by a minus sign in 
front of the third term on the r.h.s. from its counterpart considered in Eq. (5) of 

h 

Hf,u = 77^ + d^TT^ + dyTT^ + r]al3d^1T"'du'IT^ , (49) 

Mpi 

To emphasize the importance of this sign difference, we show that we recover the 
results of Ref. [H] when deriving the decoupling limit using ( 1491) . but stress that the 
Bianchi identity is then not satisfied, as a consequence of the fact that Hf^,y is then 
not a covariant tensor if ^f^^, and /i^j^ are conventionally defined. 

Up to the cubic order, the Lagrangian in the decoupling limit is then 



C = -Iw'^E-^Kp + W^X^^;} (50) 

-^ ((8ci + A)m + (8c2 - 4)[n][n^] + 8c3[n]3) + ]^h'"'xf:^ , 

with XliJ = Xlu , since both approaches only differ at quadratic order in vr, and 

x^^ = - (sci - nj, - 2(1 + C2)[u]u,, + Q - 303^ [n]\, - C2[n>^, . (51) 

Setting ci = 2c3 — | and C2 = — 3c3 + ^ to obtain the total derivative combination 
(TT5|) . we get 

x^'J = -6(c3 - \){ii% - [n]n,.) - (3c3 - 1) {[ur - [n^]) r/,. , (52) 
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which is not conserved for any choice of C3 since the reparanietrization invariance is 
not present with this choice of -ff^^, and the Bianchi identity has no reason to be 
satisfied. 

Similarly at the quartic order, we would need to impose the relation between the 
coefficients di = —6d^ — ^(24c3 — 5),d2 = Sd^ + |(6c3 — 1), ^3 = 8^5 + ^(12c3 — 1), 
and di = —6^5 — |c3, to cancel the terms of the form A^^^ddn)^. The mixing with 

the helicity-2 mode, will then enter with the quantity XfiJ as derived in [9]: 

xjS = (-i + 9c3 + 244)(nJ,-[n]ny-(9c3 + 24ci5)n^,([nV[nn(53) 
-(c3 + 8d5)vt.u {[nf - 3[n][n2] + 2[u^]) . 

As noticed in [9], not only there would then be no choice of C3 and d^ for which 
this interaction disappears, and it would always lead to higher derivative equations 
of motion, suggesting a ghost-like instability. However the fact that XfiJ is not 
conserved is an artifact of the sign choice in the expression for -ff^j^, that does not 
lead to reparanietrization invariant results. 

References 



[1] 

[2] 
[3] 

[4] 
[5] 

[6] 
[7] 

[8] 
[9] 

[10] 
[11] 



M. Fierz and W. Pauh, Proc. Roy. Soc. Lond. A 173, 211 (1939). 

P. van Nieuwenhuizen, Nucl. Phys. B 60 (1973) 478. 

H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970); V. I. Za- 
kharov, JETP Lett. 12, 312 (1970) [Pisma Zh. Eksp. Teor. Fiz. 12, 447 (1970)]. 

A. I. Vainshtein, Phys. Lett. B 39, 393 (1972). 

C DefFayet, G. R. Dvali, G. Gabadadze and A. I. Vainshtein, Phys. Rev. D 65, 
044026 (2002) [arXiv:hep-th/0106001 ]. 

D. G. Boulware and S. Deser, Phys. Rev. D 6, 3368 (1972). 

G. Gabadadze and A. Gruzinov, Phys. Rev. D 72, 124007 (2005) 



arXiv:hep-t h/03 120741 . 
N. Arkani-Hamed, H. Georgi and M. D. Schwartz, Annals Phys. 305, 96 (2003). 

P. Creminelli, A. Nicolis, M. Papucci and E. Trincherini, JHEP 0509, 003 

(2005). 

C. Deffayet and J. W. Rombouts, Phys. Rev. D 72, 044003 (2005) 
| arXiv:gr-qc/0505134] . 

L L Kogan, S. Mouslopoulos and A. Papazoglou, Phys. Lett. B 503, 173 (2001) 
tarXiv:hep-th/0011138] . 

16 



[12] M. Porrati, Pliys. Lett. B 498, 92 (2001) |arXiv:hep-th/0011152 



[13] A. Higuchi, Nucl. Phys. B 282, 397 (1987). 

[14] G. Gabadadze, Phys. Lett. B 681, 89 (2009) [ arXiv:0908llT2l [hep-th]]. 

[15] C. de Rham, Phys. Lett. B 688, 137 (2010) [ arXiv:0910.5474l [hep-th]]. 

[16] C. de Rham and G. Gabadadze. larXiv: 1006.43671 [hep-th]. 

[17] A. Nicolis, R. Rattazzi and E. Trincherini, Phys. Rev. D 79, 064036 (2009) 
[arX iv:0811.2T97l [hep-th]]. 

[18] G. R. Dvah, G. Gabadadze and M. Porrati, Phys. Lett. B 485, 208 (2000) 
|arXiv:hep-th/0005 016]; G . R. Dvah and G. Gabadadze, Phys. Rev. D 63, 
065007 (2001) |arXiv:hep-th/0008054]. 

[19] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309, 029 (2003) 
| arXiv:hep-th/0303116] . 

[20] C. Deffayet, G. Esposito-Farese and A. Vikman, Phys. Rev. D 79, 084003 (2009) 
[arXiv:0901.1314 [hep-th]], 

C. Deffayet, S. Deser and G. Esposito-Farese, Phys. Rev. D 80, 064015 (2009) 
|arXiv:0906.1967l [gr-qc] ] . 

[21] C. de Rham and A. J. ToUey, JCAP 1005, 015 (2010) [ arXiv:1003.5917l [hep- 
th]]. 

[22] A. L Vainshtein and L B. Khriplovich, Yad. Fiz., 13 (1971), 198; [Sov. J. NucL 
Phys., 13 (1971), 111. 

[23] R. L. Arnowitt, S. Deser and C. W. Misner, Phys. Rev. 116, 1322 (1959). 

[24] G. Gabadadze and L. Grisa, Phys. Lett. B 617, 124 (2005) 
|arXiv:hep-th/0412332] . 



17 



